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ABSTRACT 

The noise along axisymmetric convergent electron beams used in O-type high -power 
microwave tubes was investigated analytically as a function of several tube design pa- 
rameters. The computations and investigations were carried out for magnetically fully 
confined flow, Brouillon flow, and electrostatically focused flow for power levels near 
1 kW and perveances from 0. lx 10" 6 to l.OxKT 6 . The computed minimum noise fig- 
ures ranged between 20 and 60 dB depending on voltage, perveance, area convergence, 
and gun design. 
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NOISE IN AX I SYMMETRIC CONVERGENT ELECTRON BEAMS 
by H. G. Kosmahl and L. U. Albers 
Lewis Research Center 

SUMMARY 

The propagation of noise along axi symmetric convergent electron beams was calcu- 
lated with the density function method by assuming one -dimensional motion with a vari- 
able cross section. The linearized Boltzmann equation was solved for the alternating- 
current distribution function by assuming half -Maxwellian direct -current distribution for 
two cases: a constant electron temperature and an electron temperature proportional to 
the beam-area-compression ratio. 

The expressions for the AC noise current and electrokinetic voltage were derived 
from a Volterra-type integral equation and the noise power spectra and were then solved 
for a range of beam -gun parameters of interest in designs of medium -power and high- 
power tubes, that is, for power levels of 1 kilowatt or less, perveances between 0. 1x10 

ft 

to 1 . 0x10 , and total beam compressions up to 81. 

The complex radiofrequency -noise convection current divided by the beam compres- 
sion was shown to be invariant to area transformations, in that its application to variable- 
temperature cases leads to equations similar to those derived earlier for constant-cross- 
section, constant -temperature beams. At voltages higher than about 500 volts, the noise 
increases sharply with increasing voltage and approximately proportionately to beam 
perveances. For the same total beam compression, less final noise results if the initial 
preapertural compression is larger than the following compression. By judicious selec- 
tion of gun -beam parameters, a significant reduction of noise by a few orders of magni- 
tude can be achieved. The compiled noise figures ranged between 20 and 60 decibels, 
depending on compression ratio, perveance, and electron temperature. 


INTRODUCTION 

The nature and the propagation of thermal noise along electron beams that are used 
in O-type microwave amplifiers have been subject to numerous investigations over the 
past 15 years (refs. 1 to 7). Haus (refs. 1 and 2) has shown that the minimum noise 



figure of any linear beam -type amplifier is determined by the knowledge of two basic 
noise parameters S and n, which in turn can be determined from the noise power spec- 
tra of 1^ and Up the radiofrequency -noise convection current and the electrokinetic 
noise voltage, respectively. Analytical treatments of the noise convection current by the 
density function method, for which one -dimensional drifting beams, half -Maxwellian dc 
velocity distribution, and full uncorrelated shot noise at each velocity class (which lead to 
the problem of solving the well-known Volterra integral equation) are assumed, were 
given by Watkins (ref. 3) and Bloom and Vural (ref. 4), among others. For accelerated 
one -dimensional electron beams, Siegman, Watkins, and Hsieh, also using the density 
function method (refs. 5 and 6), carried out a numerical calculation of the noise fluctua- 
tions. To date, no' investigations on convergent accelerated multivelocity beams are 
available. The investigation of Gandhi (ref. 8), which treats the transport of disturbances 
on electron beams under magnetic compression, assumes the constancy of S and n and 
is, therefore, not valid for the propagation of noise in actual -multivelocity Maxwellian 
beams, along which both S and n vary strongly with distance. 

The advent of the space age, with its demands for ever higher frequencies combined 
with long life and the least possible power consumption, suggests the investigation of noise 
performance using convergent electron guns with axial symmetry (fig. 1). The rather 
respectable difficulties in solving a general multidimensional Boltzmann equation were 
avoided by using a "trick" frequently applied in fluid mechanics: to treat the problem as 
a one -dimensional motion but with a variable cross section. Between C and A (fig. 1) 
this approach amounts to solving the equation in spherical coordinates with the radius as 
the only variable and assuming independence of the polar angle and azimuth. The geom- 
etry of an idealized convergent gun is indeed that of a spherical cone (see fig. 1) sym- 
metric about the O-O’ axis, which goes through the center of the sphere, and with a 
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Figure 1. - Convergent electron gun. 
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solid angle r. The anode aperture A is a distance a Q away from the cathode C and, 
for a homocentric flow, all spherical shells r - x = constant form equipotential lines. 

v 

Beyond the anode aperture, which constitutes a diverging lens, the trajectories are 
no longer homocentric but change into a drifting beam converging toward its minimum 
ultimate radius b Q . Therefore, an exponentially decreasing beam radius b was assumed 
for the region between the aperture and the axial distance at which b approaches b Q . 
Smooth transition between the two regions is assured by matching the radii and their 
slopes at the aperture. 

In contrast to one -dimensional beams, along which the electron temperature T g is 
or may be assumed constant, T g increases along convergent beams. For constant cur- 
rent density across the beam, it may be shown with the help of statistical mechanics 
(refs. 9 and 10) that the transverse temperature increases proportionately to beam-area 
compression. Since * in O-type tubes - longitudinal beam fluctuations caused by the 
longitudinal electron temperature are amplified, the question arises as to the relation 
between the transverse and longitudinal temperatures. For the case of confined flow guns, 
with a strong magnetic flux cutting through the cathode, the transverse electron fluc- 
tuations are converted into longitudinal noise as a result of the adiabatic conservation of 
the electron dipole moment. Thus, in this case, both the transverse and the longitudinal 
temperatures are proportional to the beam -area compression. 

For Brouillon-type guns, or guns fully shielded from magnetic fields, little or no 
conversion of transverse into longitudinal fluctuations (in absence of collisions) will occur 
until the electrons enter the converging magnetic field region beyond the aperture. At 
this position, however, the stream has reached its final dc velocity and the velocity 
spread is small. On the other hand, all "practical” Brouillon-type guns have a small 
amount of flux cutting through the cathode so that a partial conversion of transverse into 
longitudinal noise is probable, too. 

The last case to be considered herein is the purely electrostatically focused guns 
characterized by a complete absence of any magnetic focusing fields. Again, in absence 
of true coulomb collisions, two different temperatures, T^ and T^, may exist independ- 
ently, and an increase of T t with compression will probably not lead to an equivalent 
buildup of Tj . Thus, T^ » T £ = constant seems more appropriate, and a converging 
beam with an approximately constant longitudinal temperature is dealt with. Both cases, 
that of Tj = constant and T^ proportional to area compression ratio A R , are treated 
subsequently by machine calculations and approximate closed-form solutions to permit 
extrapolations between these two temperature limits. 

The integration of the ac density function over the velocity spread leads to a Volterra- 
type integral equation for the noise convection current. The integratiqn of this Volterra 
equation with a beam having a variable density and a variable cross section requires 
knowledge of the reduced plasma frequency and, thus, the plasma frequency reduction 
factor R, The value of R has been computed for constant -velocity, constant -diameter 
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beams (ref. 11). For variable -velocity beams with changing diameters, exact expres- 
sions for R are not available, and no attempt was made to obtain an exact solution in the 
present study. However, R is not much different from 1 in the cathode -anode region, 
and the introduction of approximate values for R does not change the final result by more 
than 5 to 10 percent compared with the case R = 1. This problem is treated in appendix A. 

In order to obtain an insight into the physics of the noise phenomena in convergent 
accelerated beams, in addition to accurate computation results, the nonlinear transit 
time function in the Volterra integral equation was expanded in a Taylor series around a 
suitable center and linearized. The resulting integral equation was then solved by the 
Laplace transformation method, thus permitting a comparison with cases previously 
treated for constant-velocity, constant-temperature, and constant-diameter beams 
(appendix B). 

The computational method chosen and the problems and difficulties in numerically 
solving the Volterra integral equation were treated separately in appendix C. An excel- 
lent check for the accuracy of the computations is the amount by which the total noise 
(e.g. , radiofrequency) power S + II deviates from 1 for constant -temperature streams. 

The power series for the Langmuir -Blodgett potential is derived in appendix D. 

The significance of the investigation lies in the fact that it permits designing medium- 
and high-power O-type tubes that will have more favorable noise performance while util- 
izing convergent guns. The results of this study become particularly important for all 
cases with relatively low output powers (1 kW or less), relatively large bandwidths (FM), 
and relatively large beam compressions (high microwave frequencies), since they permit 
signal-to-noise-ratio designs of 50 decibels or more without the need of unnecessarily 
increasing the output power, even in multichannel, wide -bandwidth systems. In view of 
the fact that our ability to generate microwave power decreases with frequency as 
f -5 /2, the preceding arguments are self-explanatory. 

The International System of Units is used throughout this report. The reader is 
referred to appendix E for definitions of the most important quantities. 


ANALYSIS 

The analysis of the noise propagation problem along convergent beams with the den- 
sity function method is in the next sections and is based on the following assumptions: 

(1) The Boltzmann equation may be linearized and is collisionless. 

(2) The motion is considered one dimensional but has a variable cross section. 

(3) The longitudinal electron temperature T^ is either constant or proportional to 
the area compression ratio A^. 
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( 4 ) The direct-current potential distribution is that given by the Langmuir -Blodgett 
series. 

( 5 ) Fully uncorrelated shot noise at each velocity class is injected into the direct- 
current stream at a 90 ° angle to the cathode. 

(6) The total alternating current in the gun is zero. 

( 7 ) The velocity distribution is half Maxwellian because of an almost complete absence 
of collisions between electrons for the range of densities and temperatures applicable to 
microwave tubes. 


Variable-Temperature Case 

The analysis begins with the more general case of variable temperature. The col- 
lisionless Boltzmann equation, written for one dimension, is (ref. 6) 

df(x,u,t) _ 0f at , du 0f _ Q m 

dt at ax dt au 

where x designates the distance along the electron trajectory, 

f(x,u,t) = f Q (x,u) + f^(x,u)e^ w * are the direct -current and alternating -current distribu- 
tion functions, and the time dependence varies as e^. In the present case, since 
f Q » f p equation (1) may be linearized by neglecting products of two alternating-current 
quantities. Let E Q (x) and E^(x) be the direct-current and alternating-current fields, 
respectively. The alternating part of equation ( 1 ) will then assume the form 

9 f i 9 f 9 f i 

jct>f + u — + tjE, — + rjE — - = 0 ( 2 ) 

1 0x 1 0u ° 0u 


with 77 = e/M . 


Conventional Noise Relations 


The following conventional notations in noise terminology are used: 


Ug(a o ) = 2?jV a ; u 2 (x) = u 2 (x) + 2^= 1 - 1 


C ' — / ■ 7 — *1 ■ 0 

G?ip 2ofrpCr r p ^ a 0 a 0 


W = 0! t (u 2 - U 2 ) = 0! t |u 2 - 27jV 0 (x)j | ; U = tt £ 


1 + 2a™ w 

. V o« . 
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a r (x) 


a T =kTW = 

1 2eV„ 
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? ? i o (a o ) 

pas e u (a ) 
■os K o' 


E(x) = a - v - fe) ; ^(Xjw) * 2a T i 0 (x) • g(w) = 2a T i Q (x)S(w)e‘ 
0X 


•w 


at 


3w 


-= 2« T i o (x) [e(w> - S(w) e” w ] 


Cathode-Aperture Region 


Next, the distribution function f(x,u) is transformed into F[x, w(x,u)] using the 
method of reference 5 for the region 0 < x s a Q 


aw 

3x 3x 0w 9x 


8F, -F, 2 a 

— - 2a T r|E 0 — i - — 2. 

3x 1 ° 3w r„ 




(3a) 


and 


3u 


3F 3F 

— 2. — = 2a T u — 2. = 2G' T 2Q! T i (x) M.u 
3w 3x 3w 3w 


(3b) 
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to obtain 


jcoFj^ + u 


3F 


9x 


1 + 4“o>l (1 



i 0 (x)E 1 (x) 


§SL 

3w 


™12 
3w r 


w 


= 0 


c 1 


_x 

r 


cj 


(4) 


In a radiofrequency open circuit (total current equal to zero), the relation between the 
convection current density i^(x) and the radiofrequency field is from Maxwell's 
equations: 


R 2 (x)i 1 (x) + jcoe 0 E ^x) = 0 (5) 

An approximate expression for R(x), the plasma-frequency reduction factor, is derived 
in appendix A. From the basic definition of F^, i^ is obtained by integrating over all 
velocity classes: 


ii( x ) = -^— 
2a T 



F^(x, w) dw 


( 6 ) 


Equation (4) is a linear, partial differential equation for F^(x,w). Note that ordinary 
differential equations are obtained for Fj(x,w) (refs. 3 to 6) for constant-diameter, 
constant -temperature streams. Equation (4) is solved by writing it in the form 




3w 


-r 



ja)F 1 + 4 a 2 ?? 



-^§-E 1 (x)u(x) 
3w 1 


CO 


For compactness, 




i§- E 1 (x) = H(x,w) 
3w 


( 8 ) 


The equations of characteristics, derivable from equation (7), are 


dx 




dw 

2w 


dF x (x,w) 


P“ F .l^ + H,x.w,lr c ( l 
L u(x,w) J c \ 


(9) 
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Solving the first of the equations yields 


w =C i (1 " 


( 10 ) 


After substituting w from equation (10) into the first and third equality terms of equa- 
tion (9), an ordinary differential equation for is obtained: 


dF 


- + j — -F x + H(x,c 1 ) = 0, 


dx u(x) 


( 11 ) 


since, using equation (10) yields 


u(x,w) = (2 77V a ) 


1/2 


v» 


. + 2ct t w 


1/2 


= (2t7VJ 1/2 


V 0 (x) 


+ 2cr o c l 


1/2 


( 12 ) 


New dimensionless variables are introduced: 



o 


_ 9 _ _ _ 1 _ _ 3 _. 
dx a Q 3£ ’ 


©„ 


co • a co •' a 
o o_ 

u s (a o ) (2r?V a ) 1/2 


and the following transit time functions are defined: 



(13) 


(14a) 


(14b) 
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«(«,■") = 


(14c) 



The general solution of equation (11) may now be written with the help of equations (12) 
to (14) in the form 


F X (|,C) = c 2 e 




-f 


H(r,Cj)e 


-Je 0 [«tt,c 1 )-«(f,c 1 )] 


df 


(15) 


If c^ is replaced by its value, c^ = w/[l - ( x / r c )p from equation (10), and if equa- 
tion (15) is solved for the new integration constant c 2 , the following equation is obtained: 


c 2 


F(£,w)e 


+jO 0(|,w) +j© 0(£,w) 


+ a o e 



H(r,c)e 


■j© 0 [a(|,w)-o(?,wjj 


dr 



(16) 


where G is an arbitrary function of c^(x, w) = w/[I - (x/r c )J^. In accordance with the 
seven analytical assumptions, the initial value condition for the noise injection at the 
cathode x = 0 (potential minimum) is given by fully uncorrelated shot current in each 
velocity class w. independent of other fluctuations w. Thus, because of equations (6) 
and (16), 


F^O, w t ) ■ G(Wj) = 2a o • i 1 (0,w i )6(w - Wj) = 2a o 6(w - Wj) (17) 

where ij(0,w.) has been normalized to 1. Equation (15) is now put into the form 
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2 


F(£,w.,w) = F 1 (0,w i )e 


, ttjr 


j 3w V r 


ii(r)e L J dr (18) 

In writing equation (18), the direct-current continuity equation for i Q (x) = i 0 (a)A(a 0 )/A(x) 
and equations (5) and (8) were used. Since • Crp = ot Q a Q = 1/ 4 rjV^, equation (18) is 
put into its final form : 

f iW (i - *4 e 0 

K Sw \ r c/ 


R i(r) 


/ a O \ 2 /X ^©o[ a (S> W )-^’> W 3 

’) (1 --2-H £ 1 (£ v >e 0L J di’ 


d|’ (19) 


The defining equation is 


/ oo jq 

F^fcWj.w)® 0 


*,«) , > ie 0 s s ({) 

dw=i 1 (^w i )e 


If equation (20) is applied to equation (19), the following expression for I^(|, w j) is 
obtained: 
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2a r 


f 


ie 0 \B s U)-m,w$ V" r J Q ° 

e (w-Wj)-f- 


/ »„ \ 2 

i - - 1 


x 


1 H r ) R2(f ' ,ll<s ’ )d{ ' / 


R , x _ w i© 0 f s (lM(l,w)+3(r,w)] 

6(w)e e 0L s J dw 


/ 


~ . - w i%^ s (i)^a,wM(r,w3' 

S(w)e e L J dw) 


( 21 ) 


After the singularities are removed by applying standard relations concerning S and 6, 
the final expression for I^(|, w.) can be reduced to 

/ „ ,2 


IlU,Wj)=. 


je 0 [a s (£)-a({,wi)] x 2 


l - e. 


i<T ° / \ 2 


J R 2 (r) ^ ijdsw.) 


(■/ 


-w j®o{V^V^-[^’ w)]] 

e e Ul J d|; dw 


( 22 ) 


Equation (22) is the well-known Volterra-type integral equation of the noise current 1^ 
derived herein for the general case of a variable -cross -section, variable temperature 
electron stream with azimuthal symmetry. In addition to 1^, the second quantity that is 
needed for the calculation of the noise performance is the electrokinetie voltage U(£, w^), 
which is defined in references 1 and 2 as 


i j 0 /9 (ij) /*°° 

U(^,w i )= -e os / wFj&Wj.w) -< w>F 1 (4,w i ,w) dw 

2 “t ^0 


(23) 
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Since, for half -Maxwellian direct-current distribution, 


/ OO 

wg(w) dw = 1 


the following equation is obtained from equations (19), (23), and (24): 


( 24 ) 


U(£,w.) = -^(^w.) + 




ie 0 & s ttM(e.Wi)] 

e 



r 

t/o 


-je 0 {« s ce)-9 s (r)-[«(«.w)^(r,w5} 


we“ w dw d£' 


(25) 


From equations (14a) and (14c), the transit time functions 0 (§) and d(£,.w) contain the 
ratios V(£)/V(l) = V(x)/V a ; that is, the direct-current potential distribution in a spherical 
electron gun. This distribution has been computed by Langmuir and Blodgett (ref. 11) 
for zero velocity at the cathode. Although the velocity spread due to the Maxwell - 
Boltzmann distribution does change the potential profile in the vicinity of the emitter com- 
pared with the case of no spread, the effect is negligibly small for an anode potential 
greater than 10 volts. On the other hand, note that full account is taken of the velocity 
spread in the beam itself by the direct-current distribution function f = 2o! T i 0 (x)S(w)e“ w . 
Refer now to figure 1 and reference 13, where r c > and x = r; 

ffi4/3 (r) 

V(x) = V(r) - V(a ) (26a) 
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(26b) 


V(|) = V 

o 


•"£) 



In appendix D, o!(r/r c ) = a [j;( a 0 / r c j] is shown to be represented by the power series, 
equation (D3): 

“ r) = = " y *(* + | y| + iJ y2 ^ + °* 652 yV + ’ • (D3) 

The values for a may be obtained from table H in reference 13. 


Anode Minimum-Beam-Diameter Region 

The preceding analysis was applied to the cathode-anode region, 0 < £ < 1. Because 
the radiofrequency input of the tube is located in the vicinity of the minimum beam diam- 
eter at approximately £ = 2, both 1^(2) and U^(2) must be known to determine the 
noise performance of the tube. Beyond the anode aperture, which constitutes a diverging 
lens, the trajectories no longer remain homocentric but change into a drifting beam that 
gradually converges toward its ultimate minimum radius b Q at approximately $ = 2. 
Therefore, the following approximation for. the beam profile is made in the region 
1 < £ < 2 


bU) = b o e~ K ^‘ * (27) 

where k(£) is an appropriate polynomial in |, such that b(|) decreases monotonically to 
the asymptotic value of b Q at £ = 2. Moreover, the initial value conditions are 


b a =b ^)|=l =mb o 


(28a) 


**<*) 


1=1 



(28b) 


13 



b <%.2 =b o 

(28c) 

»’(5) 4= 2' 0 

(28d) 


In equation (28b), tan r/2 has been replaced by sin r/2. This approximation is valid 
for small angles, r < 10°. The conditions of equations (28) provide a smooth transition 
for the beam size and its slope from region £ < 1 into £ > 1. They are easily derived 
from figure 1. The condition of equation (28a), fixes the ratio m > 1 of the beam radius 
at the aperture £ = 1 and b Q at 2. From equations (28), 

b e" K ( 1 ) = mb ; - [k( 1) + «(l)lb e'*^^ = — (29a) 

r - a^ 
c o 


b o e’2 K (2) =b Q ;- [Zk( 2) + /c(2)]e" 2K(2) = 0 (29b) 

The compression ratio in the beam cross section A R , referred to the beam cross section 
at the cathode, is 


A r(£) 


b 2 (£) 



for 0 < | < 1 


(30a) 


a r U) - 


m 


b 2 (S) 


= 2 «(§>• 4 


a. 


for 1 < | < 2 


(30b) 


The transit time functions 3 g (|) and 0(£,w) from equation (14) simplify now with 
V(^) = V(l) = V a for > 1: 



(31a) 
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(31b) 


3(1, w) = 



dr 



+ 2a w • 
o 



Note that, because of the transformation of f(x, w) into F(x, w(x, u)), the direct-current 
acceleration due to E Q (which is zero for | > 1) drops out from equation (4). Therefore, 
equations (4), (22), and (25) may be used for calculating F R I R and in the region 
£ >. 1, provided that the compression ratios A R (£) (eqs. (30)) and the transit time func- 
tions occurring in them are replaced by equations (30) and (31), respectively. Thus, the 
following equation is obtained in general for 0 < £ < 2: 


j©o&s^)-^’ w i)] ( a o\ n 

I i<€* w i ) = e L Ja r^) + — - 0 o I 1 --) a r^) / 

0 \ e/ y r=0 


, 2 I i(r,w.) 

**(«> / R 2 (r) -■ --f 
1 A R (r) 


*00 

1 - / e 
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-w+j© 0 f s (^M s (r)^(rw)+a(r,w)] 


dw dr> (32) 


U x (r w i) = -I 1 (rw i ) + WjA^d) . e' 


J.©o[V*>-««^i)] 


2 , 2 
x a/. V 


1 - — 1 0 O A r (|) 


J ff o \ r c 


/ oo 

e 


r (r) 

A R ( r) 

-w+j© o [^ s (i)-o s (i')-'9(rw)+'9(r,w)] 


dw dr V (33) 


Note that, because of the character of the Volterra equation, its integration must proceed 
over the entire range 0 < | by applying the expressions for I and U within their 

respective range of validity. Equation (32) permits an important transformation: 

I (£, Wl ) is divided by A R (ij), the beam compression, and a new function J x (r w i) is 
defined as 
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(34) 


Jl&W.) s 


A r (5) 


Equation (32) can then be put into the form 


= e 




<v 2 . \ 2 /»£ 

- 1 -^) ®o f 

,a °\ r c / 4> =0 


R (I'lJxU'.Wi) 


X 



-w+j© o [o s (^)-3 s (|')-^(?,w)+5(r,wj] 

0 


d£’ dw) 


(35) 

\ 


2 F 2T 2 

H ol- ( a o / r cH 0 o ' ^ (£) is considered as a kind of a modified plasma frequency, 
equation (39) for J^(|,w.) is similar to that obtained earlier by Bloom and Vural (ref. 4) 
for the simpler case of a constant-velocity, constant-temperature, and constant-cross- 
section stream. Since the stream is still Maxwellian, all conclusions valid for 1^ apply 
as well to Jj. Especially, S + II = 1 for constant temperature beams with Ij must also 
hold for functions derived from J^. Moreover, since 1^ (in ref. 4) may be considered 
as the complex noise current divided by the unity compression ratio, is, in general, 
an invariant to compressions, but is not, at least as long as the stream is one dimen- 
sional but with a variable cross section, Maxwellian and axisymmetric 


Constant Longitudinal Temperature 

The case of constant longitudinal temperature can be easily obtained from the 
equations developed in this section. For a Q and a Q , remaining constant, equation (4) 
simplifies to 


jcoFj +u(x) 


pf 1 „ - 

-l + 4a 2 rji 0 (x)E Ax) ^ 
3x O o i g w 


= 0 


(36) 


which is now an ordinary differential equation for F^(x,u). The transit time function 
3 (£) (eq. (14a)) remains unchanged and, since V(£') s V(l) for £' > 1, equation (14c) 

S 

becomes 
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V*,w) = 



S£l + 2a 0 wY /2 

T 7 - O I 


(37) 


After integrating F^(x, w), the following equations for the noise convection current 
I, n,w.) and the electrokinetic voltage U, (£, w.) along a beam with a constant longi- 

VJ X X « V/ X 

tuainal temperature and a variable cross section are obtained: 


I l, 0 (5.V = eie °^ S<{Ms<l ’ W ^ + ^ (»—) ®o f l R 2 (S’)l 1;0 (i’,w l )- 

OX C / o 


A r(4’) 


r 

poo 

Vi 


-w+j© 0 [a s (|)-^ s (f )- 1 9 o (|, w)+a o (|% wj] 


dw d£’> (38) 


, x , . JdA(S)-a o « f Wif| X 2 f a \ 2 r 

U l,o<«’ w i> = - I l,o(^ w i> +w i e J - ^ (j “j 0 o J 


r |'=0 


R 2 (|') 


XI l.o^’’ w o) ‘ A R ( ^ ,) # we " e 

'0 


«/■ 

Jo 


-w. 


i0 o [? s (s)“V*')“V& w ) +5 0 ^’’ w 3 


dw d^ (39) 


Noise Power Spectra 

At this point in the analysis, the noise power spectra can be calculated by averaging 
|U| 2 , |l| 2 , and I - U* over the half -Maxwellian distribution (ref. 1). From reference 1 
are obtained 


^(|) = <|l(^w i )| 2 > 


|l(|,w.)| 2 e *dw. 


(40) 


<5 - W . 

$(l) = <|U(£,w.)| Z > * / |U(4,w.)| 2 e x dw. 

'0 


(41) 
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(42) 


S*O0 

n(!) + iM& = (1(1, Wj> • u^w.)) S J I(?,w i )U*(|,w i )e’ Wi dWj 


s«) -[*«>• $a)-A 2 (^)] 


1/2 


(43) 


In equations (40) to (43), either 1^ and or I 1 Q and should be used. Haus 

(refs. 1 and 2) has shown that the minimum possible noise figure of any linear electron 
beam amplifier is given by the expression 


T 

F min<^ = 1 + [ S (£> - — ( 44a ) 

T o 

Of particular importance to the present analysis is the value of F m i n (l)^ = 2 > which, 
after it is multiplyed by kT Q times the bandwidth, represents the least possible noise 
power at the signal input of the amplifier. The cathode temperature T is equal to the 
physical temperature of the cathode surface. Note that, in the present definition, both 
S and II are normalized dimensionless noise spectral functions. 

Equation (44a) is valid strictly for lossless circuits. For practical interaction cir- 
cuits, (S - 33) in equation (44a) should be multiplied by the factor |l + (d/xjJ, where d 
and Xj are the loss parameter and growing wave parameter, respectively. For narrow- 
bandwidth (e.g., < 5 percent) high -impedance -interaction circuits, such as cavities, 
this factor is very close to 1. For wide-band tubes (e.g. , helical circuits), the factor 
may be as high as 1. 2 to 1. 5 in the frequency range of interest (< 12 GHz). 

For F m - n » 1, as is the case in power tubes, a very accurate approximation is 
obtained by neglecting 1 in equation (44a). Thus, 

“ [ S <«> - n «>] ^ < 44b > 

RESULTS AND DISCUSSION 

Choice of Parameters for Medium- and High-Power Tubes 

Before the computational and analytical results are reviewed, the meaning and the 
range of beam -gun design parameters are discussed. It should be borne in mind that the 
selection of these parameters is dictated by power tube requirements, such as efficiency, 
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gain, and beam transmission, and not by minimum noise requirements. Therefore, the 
results of this study, if carried out over a suitable range of typical numbers, show 
representative noise levels associated with power tube designs. The beam -gun param- 
eters used in this study are © . a /r , x . cr, b , and m. 

Cr o c*. a* ’ o’ 

By definition, 0 Q - a> • a Q / (2?jV 0 )V2 = (a)/u Q )a o is the transit time in radians (in 
the absence of space charge) of an electron having the full anode velocity u Q - 2rjV^^ 
all the way between cathode and anode. The actual transit time is, of course, larger 
than © Q and is given by O^w^ • © o< However, 0 Q is proportional to the distance 
cathode -anode a Q and is therefore an important geometric gun parameter. Guns in 
medium- to high-power tubes are designed, in general, for 0 Q w 4 to 12 radians. The 
approximate range of a /r , the distance cathode-anode divided by the radius of curva- 
hire of the cathode, is 1/3 to 2/3. 

The normalized plasma frequency at the anode aperture is X = a) _/&>. The prac- 

a pj a 

tical range for x a is 0. 01 to about 0. 2, which corresponds approximately to micro - 
perveances of 0. 1 to 1. 0 when the final beam size b Q is referred to. Notice that Q , 
the plasma frequency in the interaction region, is larger than a) by the factor ’ 

m=b a /b o- 

In medium- to high-power tubes, Q /ci) « 0. 1, and the approximate relation holds 
for (u>/u Q )b 0 = 1/2: 


R o X o = W loo x pj/ 2 ; x 0 = ^-2 = mx a « 3. 5xl0 2 x P^ 2 (45) 

The quantity cr = kT /2eV = kT /2eV . If T is chosen to be 1100 K, which is close 
to the true cathode temperatures, the anode voltage V Q = V a may be expressed by a Q . 
Since, in tubes designed for about 1 -kilowatt power, V Q = 1000 to 10 000 volts, the cor- 
responding a Q range is 10”^ to 0. 5x10“^. 

The ultimate beam radius is b Q , and the beam radius at the anode aperture is b R . 
From equation (28a), b a =b(^)^_^ = mb Q . In modern, high -efficiency tubes, 

(or/u o )b o « 0. 5; therefore, 


b « 
o 


u u 

0. 5 - 2 . = 0. 5 — 2 — 
u) or • a Q 




The distance a Q itself may be determined from the relation 




Cl) 


1. 29X10 5 



(46) 
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where a Q is in meters, © Q is in radians, and oo is the angular frequency. 

The total area compression ratio A R is equal to m 2 /^ - (a Q /r c y| 2 . In this study, 
a Q /r c = 1/2 and 2/3, and m = 2 and 3 were used to yield an A R from 16 to 81. Notice 
that a given A R can be achieved by making the compression y [i - (a Q /r c )] 2 in the 
spherical diode larger (smaller) and the subsequent compression m 2 smaller (Larger). 

Two frequencies, 2 and 10 gigahertz, covering the range of interest to space com- 
munications and broadcast satellites, were selected in this study. 

Some of the selected beam -gun designs are certainly unrealistic. If, for example, 

© Q is relatively small and the frequency is high, according to equations (46) and (47), 
both b Q and a Q would be small. Because monotonic beam convergence is required with 
boundary values specified by equations (28), an unrealistically rapid convergence is 
required immediately following the anode aperture, that is, at 1. 0 £ <, 1, 2. Obviously, 

a gun designer would not select such an odd combination of parameters; this compilation 
was made only for the sake of completeness. 

The relations (eq. (45)) between the actual tube perveance P , the anodal compres- 

2 o ^ 

sion factor m , and the relative plasma frequency x a are summarized in table I. In 
deriving the numerical values, (o>/u 0 )b Q = 1/2 was assumed. 


TABLE I. - RELATIVE PLASMA FREQUENCY AS 
FUNCTION OF BEAM PERVEANCE AND 


ANODAL BEAM COMPRESSION 


Actual tube 
perveance, 

P M 

Anodal 

compression 

2 

factor, m 

Relative plasma 
frequency, 

* 

0. IxlO' 6 

4,9 

3.1X10 -3 , 1.36X10 -3 

1/4 

4 

3/4X10" 2 

1/2 

4 

3/2 

1 

4 

3 

1/4 

9 

1/3 

1/2 

9 

2/3 

1 

9 

4/3 


Discussion of Results 

The results of this study, plotted in figures 2 to 14, were obtained by solving the 
expressions for I, U, S, and II on a high-speed, large-storage IBM 7094 computer. 
Figures 2 to 5 show the quantity S - II plotted against £, the normalized distance, 
between the cathode (| = 0) and the anode (£ = 1) and at £ = 2, where the compression 
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is complete for the variable -temperature case. The graphs cover the following range of 
parameters: a Q = 10“ 3 , 10“ 4 , and 10“®; A 3 = 10“ 3 , 10“ 2 , and 3xl0“ 2 ; 0 Q = 4 and 8; 
a Q /r c = 1/2 and 2/3; m = 2 and 3; and oj = 2 and 10 gigahertz. Examination of figures 2 
to 5 shows an approximately exponential rise in noise with distance between P = 0 and 

O A 

1 = 1. Beyond £ = 1, the low-voltage cases (a Q = 10“ and 10“ ) exhibit a leveling off of 

S - n, which indicates that S and n are both approaching their respective asymptotic 

values at large distances from the cathode. In contrast to this behavior, the high-voltage 

cases (a = 10"®) continue to rise steeply with distance beyond | = 1 for high perveance 
u o -2 o o 

values (A & s 10 ) and less steeply for the low perveance value (A a = 10 ). Notice that, 

at | = 2, the values of S * II in high-voltage cases are several orders of magnitude 

higher than those in low -voltage cases and that, for the same cr value, an order of 

magnitude variation in S - II may be found depending on the size of A: and a /r . The 

curves show a similar behavior for 0 < 4 < 1 and branch off for different values of m 

(e.g. , different postanodal compressions). Also, the exponential rise in noise for 

0 < 4 < 1 is in agreement with analytic predictions derived from approximate closed - 

form solutions in appendix B . 

The results of figures 2 to 5 are summarized in figures 6 to 9, where S - II at 
4 = 2 (which is of interest to tube designers) is plotted against 1 /ct q (i. e. , against a 
quantity proportional to voltage, with m, A„, 0 , and a /r as parameters). An in- 
crease in noise with increasing voltage is observed in the range l/a Q > 10 (i. e. , for 

V > 500 V), a range of practical importance to designers of tubes with power levels of 

W Q 

about 1 kilowatt. (Note that the value l/a Q = 10 corresponds to a value of V Q of 

approximately 50 V and is therefore only of academic interest. ) The larger the values 
2 

of A a and A R’ the larger is this increase in noise. A significant reduction in noise is 
possible at l/a Q = 10®, that is, V Q « 5000 volts and above for A 2 10“ 3 . From table I, 
the latter is seen to correspond to values of P < 0. 1x10“®, which are smaller than con- 
ventional values. Another significant conclusion can be drawn from a comparison of 
figures 7 and 8. The results summarized in both figures were compiled for the same 
total compression ratio of 36 but for different preanodal and postanodal compressions: 
in figure 7, the initial compression in the spherical part of the gun is 4 and the final com- 
pression is 9; in figure 8, the sequence is reversed. For higher voltages (e.g., 

^o > 10®), the resulting (S - II)^ is seen to be substantially reduced if the initial com- 
pression a /r is made larger and the subsequent compression smaller. Physically 
speaking, the "cooling" of beam noise is done more effectively in the low -voltage, large- 
velocity -spread region than could be achieved if the noise, due to rising temperature, 
were injected in the high-velocity, negligible -spread region Thus, guns with perveances 
> 1/ 4x10“ should be designed for high initial compression to be more effective in 
reducing noise. 

The preceding discussion is valid, strictly speaking, for fully confined-flow beams. 

As already mentioned in the INTRODUCTION, less noise should be expected for Brouillon 
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compression 

factor, 



Normalized distance, £ 

(a) Relative plasma frequency, xf, 0. 001. (b) Relative plasma frequency, 0. 01. (c) Relative plasma frequency, xf, 0. 03. 

Figure 4 - Quantity S -IT as function of normalized distance for variable-temperature case. Transit time parameter, 0 O , 4; ratio of cathode-anode distance to cathode radius, a 0 /r c , 1/2; noisiness, 
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Figure 8. - Quantity S - IT at normalized distance of 2 as function of l/o (pro- 
portional to voltage) with various values of relative plasma frequency (pro- 
portional to perveance) for variable -temperature case. Transit time param- 
eter, 0 O , 8; ratio of cathode-anode distance to cathode radius, a Q /r c> 2/3; 
anodal compression factor, m, 2; area compression ratio, A R , 36. 



Figure 9. - Quantity S - 17 at normalized distance of 2 as function of l/o (pro- 
portional to voltage) with various values of relative plasma frequency (propor- 
tional to perveance) for variable-temperature case. Transit time parameter, 
0 O , 8; ratio of cathode-anode distance to cathode radius, a Q fr c , 2/3; anodal 
compression factor, m, 3; area compression ratio, A R , 81. 
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focused beams. Although no detailed investigation of this case was made in the present 
study, it is felt that the difference would not be as large as could be expected for two 
reasons: First, in the absence of magnetic fields at the cathode surface, some electrons 
will leave the surface not in normal directions, but rather at all angles between 0° and 
90°. This angular spread, which is almost entirely inhibited in confined-flow guns, was 
not included in the present study, although it certainly contributes to additional noise. 
Second, an entry into a strong magnetic field always introduces some noise in the high- 
velocity region where cooling is ineffective. 

The discussion now proceeds to figures 10 to 14, compiled for constant-temperature 
(CT) "ideal" electron beams, ideal meaning the assumption of a flow in which all elec- 
trons leave the cathode normal to its surface. Again, this assumption is justified only 
for confined-flow beams, and therefore in purely electrostatically focused beams (along 
which constant temperature may be expected only for this ideal situation) will have higher 
noise in practical cases than that plotted in figures 10 to 14. However, the calculated 
numbers still represent a useful lower limit of noise. A review of the plots for constant - 
temperature cases reveals the expected result that the (S - n)^ numbers are two or 
three orders of magnitude lower than those computed for variable -temperature cases at 
higher £ values. At lower £ values (£ < 1/2), the difference is small because not 
much noise conversion had taken place yet. Qualitatively speaking, the general para- 
metric behavior is not different from that discussed for variable -temperature cases: 
figures 10 to 14 show the beneficial effect of higher initial compression (a / r large) as 
well as the detrimental effect of higher perveance (X„ large) on the magnitude of (S - II). 
Also, a general increase in noise with increasing voltage is clearly evident. 

The discussion of the beam -gun parameters is now completed with an examination of 
the effect of the gun length © Q on (S - n)^. A comparison of figure 4 with figures 3 
and 5 shows that the effect of © Q is not very profound and is only of the order of 20 per- 
cent or less. Longer guns (0 Q ^ 8) seem to be better for high -perveance designs, and 


Reiative 

plasma 



s-n 



Anodal compression factor, m 2 


6 1 I l l l I I I l I 1 

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 

Normalized distance, 8; 

Figure 10. - Quantity S -17 as function of normalized distance for constant-temperature 
case. Transit time parameter, 0 Q , 8; noisiness, a, 10~ 3 ; ratio of cathode-anode distance 
to cathode radius, a Q /r c , 1/2. 
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Anodal 

compression 



Figure 13. - Quantity S - TT as function of normalized distance for con- 
stant-temperature case. Transit time parameter 0 Q , 4; noisiness, a, 
1(T 5 ; ratio of cathode-anode distance to cathode radius, a 0 /r c , 1/2. 


shorter guns are slightly better for low-perveance designs. A similar behavior may be 
deduced from constant-temperature cases by comparing figures 12 and 14 (for © Q = 8) 
with figures 11 and 13 (for 0 Q = 4). The behavior of either unrealistically short guns 
(0 Q < 2) or very long guns (0 Q > 12) was not investigated in detail. However, a few com- 
putations made for 0 Q = 2 and 12 indicated a sharp increase in noise for 0 Q = 2 com- 
pared with cases for © o = 4 and 8. Additional noise not considered in this study may also 
be caused by ripple or scallops on the beam (particularly on beams other than confined 
flow) as a result of improper launching or transition through electric or magnetic lenses. 

To conclude the discussion, the error resulting from the omission of the cathode - 
potential minimum region in the Langmuir -Blodgett equation will be estimated. This 
omission is equivalent to the assumption 5 (p. 5) which lets the potential minimum co- 
incide with the cathode itself and thus sets the injected noise current equal to the full shot 
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Normalized distance, E, 


Figure 14 - Quantity S - TT as function of normalized distance for constant-temperature case. Transit time 
parameter, 0 Q , 8 ; noisiness, o, 10 r. 


noise. From the exact analysis of the noise behavior in the vicinity of a potential mini- 
mum (ref. 12), it is known that the shot noise smoothing factor r is approximately 0. 3 
at radiofrequencies f smaller than the plasma frequency f m at the potential minimum 
and approximately 1 at f >fp m - Since in high-power microwave tubes, f > 3 giga- 
hertz, the results of this study may be too high by 5 decibels or less at microwave fre- 
quencies below 3 gigahertz and quite accurate above this value. Moreover, it is the 
range of very high frequencies (e.g., the 10 to 35 gHz range) for which the present study 
is of special importance because of difficulties associated with generation of power and 
high-compression ratios. Note, also, that substantial smoothing of noise takes place in 
the low velocity region following the shot current injection and that the magnitude of the 
potential minimum being less than 0. 5 volt is negligibly small in comparison with the 
kilovolt anode potential. 
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SUMMARY OF RESULTS 


The following results were obtained from the study of noise propagation along axi- 
symmetric accelerated electron streams for two cases: a longitudinal electron temper- 
ature proportional to the beam area compression ratio (confined flow) and a longitudinal 
electron temperature constant throughout the entire compression region (ideal Brouillon 
and electrostatic focusing), in both of which flows all electrons leave the cathode nor- 
mal to its surface: 

1. For voltages above a few hundred volts, the noise rises sharply with increasing 
voltage and less steeply with increasing perveance. 

2. For the same total compression ratio, stronger initial compression in the 
spherical part of the gun followed by smaller subsequent compression produces substan- 
tially less noise than the reverse situation, especially at larger perveances, 

> 1/4X10" 6 . 

3. For conventional gun designs, with a transit time parameter 0 Q = 4 to 10 radians, 
the noise figure changes only very weakly with © • however, the noise rises sharply at 

R O 

higher perveances, > 1/2x10 , with decreasing gun length, 0 Q < 2. 

4. The introduction of a new variable, the complex noise convection current divided 
by the local compression ratio, showed that the basic equations for noise current and 
electrokinetic voltage become invariant to area transformations, thus reducing the prob- 
lem to that of one -dimensional Maxwellian streams. 

5. The computed noise figures are expected to be fairly accurate for the case of 
confined-flow streams, the most important being high-power, long-life tubes, and are 
expected to represent a lower noise limit for Brouillon flow and electrostatically focused 
beams. 

6. For the range of parameters applicable to tubes designed for power levels of 
approximately a few kilowatts, the computed noise figures ranged between 30 to 60 deci- 
bels in variable -temperature cases and between 20 to 30 decibels in constant-temperature 
cases. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 5, 1969, 

160-21-06-02-22. 
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APPENDIX A 


APPROXIMATE EXPRESSION FOR PLASMA-FREQUENCY-REDUCTION FACTOR 

In dealing with electron beams with finite cross sections surrounded by a conducting 
boundary such as a metallic tube of radius a, the effective plasma frequency is 

= R • o>p, where R is the well-known reduction factor. For beams having a con- 
stant radius b and moving with a constant velocity u Q , R = Rj(w/u 0 )b, (b/a) has been 
computed and plotted for various values of b/a and (to/u Q ) • b(e.g., see ref. 11). 

Notice that R approaches 1 for very slow beams, u Q small, or for very high frequencies 
to; and R — 0 for the reverse conditions. 

The reduction factor R has not yet been computed for beams with a variable velocity 

u Q and/or a variable radius, b. As shown in figure 1, both beam parameters u and b 

change simultaneously. However, an approximate solution consistent with the assumption 
of a homocentric motion with variable cross section may be obtained: Write r for x 
and assume (see INTRODUCTION) independence of the polar coordinates cp and Q for i 
and E (i . e . , i = i(r) and E - E(r), only). The continuity equation is 

_L JL(r 2 i) = (Al) 

r 2 0r 0t 

For the Poisson equation, write (ref . 14) 

JL _L (r 2 E) -R 2 -£- (A2) 

r 2 5r 


where R = R(o>b/u) but not of the coordinate r itself, although both b = b(r) and 
u = u(r) depend on the position. The formulation equation (A2) is standard for cylin- 
drical beams of constant radius. The integration of the direct-current case, 3/ at = 0 in 
equation (Al), results in 



which is identical with the direct-current continuity equation. Partial differentiation of 
equation (A2) with respect to t gives 
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or 


= = (r 2 i) 

r 2 9r \ at/ e 0 at ~ e Q r 2 0r 


®M-j„E = -5!i(r) 

3t 


(A3) 


Then, from the earlier notation r = x (eq. (5)), 

R 2 pJ*M|i ( x ) + jco€ 0 E 1 (x) = 0 
j_u(x)J x 

Equations (A3) and (5) are believed to be fairly accurate as long as the changes caused 
by the fluctuations of the radiofrequency field at any instant are larger in effect than 
those produced by the ordered electron motion. From reference 14, the approximate 
expression for R is 


R 


2 




(A4) 


Now, from figure 1, assuming that the shield radius a « r at £ « 0 and b/a « 1/2 
for ij 5: 1 yields 


b(x) = all 


i = a 




b(a o )=a^-!°j (A 5) 

Note that a, the shield radius, must not be confused with a Q , the distance between 
cathode and anode (see fig. 1). Substituting equation (A 5) into equation (A4) gives (for 
0 ^ {s 1) 
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V(€) 

V(a Q ) 
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1 + 


3 (i-V 


a 


4/3 


f 




4 b \> a v 3 (S 


i + l/!aA + Ap>{' 

2 Vc ) 2 Vo , 


l/ a o 


l/^o 

2 l r 


a 3 + . 


n2 


(A6) 


Notice that R(0) - 1 and 0 < R(l) < 1, and 0 a = oj/u(a o ). 

The usual design of microwave tubes in the interaction region is such that 
a)b 0 / u 0 « 0. 5. Since the beam radius b Q at the anode aperture is about two to three 
times b Q , the value in the interaction region, and since u a = u Q in most cases, 
cob a /u o will be approximately 1. 5 to 3. 0, and R from expression (A6) assumes the 
value R(l) « 0. 7 for £ = 1. 0. 

At low values of £ < 1/4, where all significant changes in equations for I(£, w.) 
and U(|,w.) and S, II, and A take place, V(l/2) is «0. 25V a and R(l/4)«0. 94. 
Thus, R is close to 1 where a more accurate knowledge of it is required. 
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APPENDIX B 


APPROXIMATE SOLUTIONS BY LAPLACE TRANSFORMATION 
Longitudinal Temperature Proportional to Beam Compression 

Because of the nonintegrability of the transit time functions 3(|, w) and their non- 
linear characters and the presence of the terms R^(|), |i - (a Q /r c )^j, and 
|I - ( a 0 / r c )fj i n the kernel, the Volterra integral equation (22) cannot be solved by either 
one of the two methods most commonly applicable to this type linear integral equations, 
for example, expansion of the kernel in eigenfunctions, if they exist, or use of the Laplace 
transformation. The latter method was successfully employed by Bloom and Vural 
(ref. 4) in solving the noise propagated on a one-dimensional, drifting beam. They ap- 
proximated their transit time function for low values of a and 0 Q as a linear function 
of distance for a drifting beam: 



exp(jow0 Q ) 


(Bl) 


In an attempt to gain a qualitative insight into the physics of the phenomena associ- 
ated with noise propagation on accelerated variable -density beams, an approximate 
analytical solution will be developed and exact computations made, with the following 
simplifying assumptions: 

(1) The beam is only slightly convergent (e.g. , a 0 < r c ) so that y < 1 and 
a(y£) « -y £ jl + (4/5)y^j according to equation (D3). 

(2) The potential profile is only slightly nonlinear. 

(3) The measure of the relative velocity spread a Q is small compared with 1. 
Equation (14b) is approximated by 



(B2) 

Because of assumptions (1) and (2), the integrand in equation (B2) is expanded around a 
suitable expansion center in a Taylor series. Thus, 
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H yS ) 




(B3) 


with 


f l<U = 






2/3 


W = - 


5 r 


3 . , a v 5/ 3 


(B4) 


Equation (B2) is integrated to the first two terms and yields 


m = 



2/3 


w- - 


(B5) 


Similarly, the integrand of equation (14c) is approximated and expanded into a Taylor 
series (y = a Q /r c ): 



(B6) 


(B7) 
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w = -ay 


a o wa f 3 * | ft - y^p ) 3 fr ^) 1/3 (* + f y ^°) 1/3 (* * ~ y ^o) 


y«o ( ! 


l+*jd* /3 (1- JU 2 + 2a?*J* /2 


The approximated expression for the transit time function is then 


3(1, w) = «| /3 f 2 (g- s + W^-v « 


The expressions (B6) and B9) are linearized by setting £ = | Q * £. Furthermore, for 
ct q « 1, because of assumption (3), the radicands in equations (B7) and (B9) are expanded 
by taking the first term only. In spite of the fact that 0 < w < °°, equation (B9) is a good 
approximation because of the presence of the weighting factor e~ w in equation (22). 

With the aforementioned simplifications, the following equations are obtained for the 
difference of the two transit time functions: 


V«> - w > = 2 ~ v* l 1 - ■ «<> - - 

y hi v 15 


= 


(BIO) 


2 — (it— y| 0 - + 

y2| 2V 15 


(BH) 


Similarly, for the difference of 3 g (£) - 3(|, w), resulting from equation (31a) minus 
(31b), the following equation is obtained for the range 0 < £ < 2: 


38 



°o w 


m 




e +K (4 , )l' dr * CT w __HL! 





dr 





_jl_ / e K( y^ 
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(B12) 


with the choice of a suitable 1 < £ Q < 2. 

Now, a new dependent variable is introduced, J(£, Wj) = I(|,w.) times the beam com- 
pression (A r ) _1 defined in equation (34). Thus, from equations (22) and (32), 


I l(^ w i ) (l - — I ] =J 1 (^ w i ) 


(B13a) 


liU.Wj) • e- 2 ' «• *«> (l - !H.\ A = J 2 ( { , Wl ) 

a V m 


Equations (22) and (32) can then be written in the form 


r^Wj) = e^ ia ° Wi ^ 0 ° + hL ^1 - 0 o R 2 (^ o ) y Jia% w i) 


X 


[-/ 


-w j 0 o^l a( ^" r)w . 
e e dw 


(B13b) 


dr (B 14a) 
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m 


j 2 ^» w i) = e 


2 a o w i 


0 


o 




• r 2 <«o> 





X 


1 - 




-w+jejjjo- (|-4')w 

e 010 dw 


+ ®y<t 0 > 



-W+j© 0 ^ 2 <T 0 (|-|’)w 

e 


dw 



j 2 (i*,w.)dr 


(B14b) 
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Notice that R (£) was approximated under the integral by suitable averages of R (£ Q ), 

since R(£) is a slowly varying function over £. The Volterra integral equations (B14a) 

and (B 14b) for the new function J 1 (|, w^ are similar, as already mentioned, to the 

Volterra equation obtained earlier by Bloom and Vural (ref. 4) for the simple case of a 

constant-diameter, constant-velocity, and constant -temperature beam. Note that, in 

the present case, the exponentials are modified by the appearance of j8 > 1 and that the 

2 

relative plasma frequency A a is also modified by the factor 



0 o r 2 <«o> 


Equations (B14a) and (B14b) can be readily solved by the Laplace transformation method, 
for which is written 


J(4,w i ) = e 


j0W 0 o ,72 


+ \X 


a 






(B15) 
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a 


(B16) 
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m - ?) = +jo 0 ' 3 


LzS. 

i - wo 0 u - r) 


%L /o 


e e aw 


(B17) 


Application of the Laplace transformation and of the convolution integral theorem to 
yields 




F(p)=- — x|F(p)- 0, 


P - l^oOoWi 


T_je-L 
01 1 


d£ 


r (bis) 


F(p) = 


(P - j0i CT o 0 o w i) 


A 1- /a 

i + i40 n — -i e 


j(p//3 !©(,%) Z* 00 


t 2 u \ p 0,0 


W Z* 00 e -v 

O «/p/^i CT 0 e 0 V 


dv 


_o 

Except for the presence of the factors /3^0 q and X a , the algebraic equation (B18) for 
F(p) is identical with that obtained and discussed earlier by Bloom and Vural for the case 
of a nonconvergent, drifting beam with constant electron temperature. The quantity 
F(p) has two purely imaginary poles 


P = P 4 = 


lV"i 6 0 


and 


P = P s = -3S^ 1 cr o 0 o 


where s is the root of the equation 


E-4316 
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The nature of the solutions and their interpretations are described in reference 4; 

therefore, the results are only summarized herein. The two poles represent the two 

simple propagating waves with amplitudes A and B, which are given by the residues. In 

the present case, interest is in medium -power tubes, with p ranging from a few hundred 

£* 

to a few thousand watts, and in conventional perveances Rn < 10" . Consequently, X 
will be approximately 0. 1, a will be between 10"^ and 10“°, and x/a » 1. For the 
condition x/o » 1, the amplitude A of the kinematic wave with the propagation con- 
stant p. = ji3 1 e o a o w i dies out to zero, and only the space -charge wave B survives 
with distance. In addition to the two propagating waves, there is also an aggregate of 
nonpropagating, Landau -damped waves. The existence of this aggregate may be con- 
cluded from the requirement that the vector sum of the amplitudes of all waves at the 
origin must be 1 so as to represent the unity -injected (normalized radiofrequency) cur- 
rent. Let the aggregate be represented by 


F(l-) = 


F e 
n 




n 


F(0) = 1 
F(«o) = 0 


(B20) 


Since the discussion in this appendix is concerned with the approximate analytical 
behavior of the noise problem, a sufficient description is obtained with one nonpropagating 
fast wave only : 


F({) a F /° 0 ° SOr l'“l > = e °o®o*O r l-“l> 


(B21 


The amplitude of the slow propagating wave is 


B = ResF (p) =lM. 

P = P S h'(p) 


(B22a) 


P=Pc 


With 
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and 


one obtains 


g(p) 


P=P S = - l^lW 1 - [-1 Wl (s + w l>]" 1 


h'(p) =*M 

u/ n=n 

v H S dp 


P=P C 


x 2 L jp/^i e„a ( 

h(p) = 1 + j — 0 O (— - j - 

•2 V * 


3.0 a /»<*> \ 

±11 f £!dv) 

l 0 o «^p//3i© 0 ct 0 V J 


(B23) 


(B24) 


B = 




(s + w.)l-| - /3 1 s 2 ’ 

Vo 


(B22b) 


From equation (B19), s is found from integration by parts of f°° (e~ v /v) dv to be 

— (— ^ < B25 > 

h Vo / 


The expression for J^(£, can be written in the form 


-js/3.,0 a Q | CT n 0 o^ r l _O! l^ 

J 1 (?,w.) = Be 10 0 + (1 - B)e 0 0 1 1 


(B26) 


For the desired noise current I^(£, w^), the following is obtained from equation (B13a): 

•WAV . i - b “o e o° r r“i ){ 




B 




K*) 


for 0 < 4 < 1 (B27a) 
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IjU.Wj) =- -e 


mi_ e -[ 2 /< r c- a o>l *«({)• i Be - 1S V, 


a 


i- 


o 


a0(jr 1 -o! 1 )| 

+ (1 - B)e 1 1 for 1 < 4 < 2 (B27b) 

The boundary conditions for I^(o, and 91^(1, w^)/d£ are 1 and 2(a Q /r c ) 

+ 5/30 o %Wi, respectively, as may be verified from equation (22). The following equation 
is obtained from equation (B27a): 


k 1 (|,w i rj 

. _L = 


2yB - SfVjV + (1 - B)a o e o (jrj - a,) 


2 — + )' 3 l e o°o w i 
r o 


Bs + w. 

r l = h 


a l = ~ 


1 - B 


2y 

cr 0 
o o 


(B28) 


For A / a » 1 applying equation (B25) to equation (B22b) gives 
a o 


B 


iL W 1 L <Wo 

~~A * ) A x a Av 


(B29) 


Therefore, r^ « ^(s + 2w i )0 Q and 


, ,, „ , B .-WlW . 1 - B 

= e + e e 


(i - y£) 


(1 - y£) 


(B30) 
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If S * (l/^X J Cq is taken for X a /a Q » 1 and substituted into equation (B30), one 
obtains, for the more important range of < 10, 




e - 3 V^l x a 0 o* + e 2 ^ 

2(1 - y£) 2 2(1 - y£) 2 


e y£ 

(iV^W y ) 5 ^ -(iV<WV y )« 

“ T 6 

(i - yl> 2 

2 


(1 - yl) 2 


cos 


(#l x a ®0 - iy )« 


- — [cos coshyl + j sin sinhyjl (B31) 

(i - y«r L J 

Equation (B31) is now compared with a corresponding case 1^ Q (|, w^, obtained from 
equation (38) for a drifting, nonconvergent, constant-temperature stream (/3 = 1, y = 0, 

T g = constant): 


Ii Q (^w i ) -« cos y| + j -w t sin y| 

i,o i x 1 

w cos y^ for -« 1 (B32) 

X 

Thus, a rising electron temperature produces complex current waves whose amplitudes 
grow even stronger than the beam compression ratio, that is, stronger than the distance 
squared, and are complex even if the contribution of Wj terms in the exponent and in B 
of equation (B30) is neglected. The factor -^3 > 1 reflects the effect of an accelerated, 
but nonconvergent, flow. The effective plasma frequency appears larger by a factor 
yffi than that in drifting beams. This result may be understood from the nature of the 
Lagmuir -Blodgett potential shape, which rises slowly over a relatively large distance. 

A corresponding similar expression can be obtained for I^(£, w^) in the range $>1 
from equation (B27b). Also, by substituting the equations for Ij(|, into the expres- 
sions (25) and (33) and by approximating the transit time functions 3(4, w) in a manner 
developed in this appendix, an approximate equation for U(£, Wj) could be obtained and, 
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from the knowledge of I and U, S and II, could be calculated. However, since U is 

/ 4 

proportional to powers of (A/ a) , the accuracy of the approximations would no longer be 
satisfactory. In contrast to U » 1, the magnitude of J-^(£, w^ is of the order of 1, and 
the approximate expressions (B30) and (B31) are fairly accurate. In spite of the difficulty 
of producing a simple and accurate expression for U, the general, qualitative behavior 
of noise can be predicted to rise with distance much more sharply than constant- 
temperature beams. This prediction follows from the fact that U increases with 1^, 
and 1^ grows stronger than the distance squared. Thus, in agreement with experi- 
mental findings, much higher noise figures should be expected along magnetically com- 
pressed beams. 


Longitudinal Temperature Constant 


The second case of interest is the converging beam, moving through same potential 

profile as in equation (Bl), but having a constant longitudinal temperature. The analysis 

of such a case would reveal the effects on the noise behavior of beam convergence and of 

the direct-current potential alone. The analysis begins with the expansion of the transit 

time functions 3 (£) and 0 (£,w) (eqs. (14a) and (14c)) in a similar fashion as done in 
s o 

equation (Bl). A transit time parameter /3 is obtained: 


0 O = 2 


a 


2 r 


.2 ,2 




(B33) 


With the help of equations (B16) and (B17), equation (38) may then be written as 


1. e 


(?,w i ) = 


+ ir 2 

3 + ] A„ 




.(r^) • m - r) 
(i - y£') 2 


d£’ 


(B34) 


2 

with K(| - $’) given by equation (B17) and |3 = /3 Q . Again, R (£’) has been approximated 
by a suitable average R 2 (4 Q ). In its form, equation (B34) cannot be solved by the Laplace 
transformation method. The solution is achieved by expanding (for y£ < 1) 


- = 1 + 2y£ + 3y 2 £ 2 + 4y 3 £ 3 + . . . 

(1 - yl) 2 


(B35) 
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Dropping temporarily the terms containing £ in equation (B35) makes equation (B33) 
solvable by Laplace transformation. In analogy to the result obtained in equation (22) for 
\/a » 1 is 


V ( 5 ’ W i ) = Be 




+ (1 - B)e 


+ja o 0 o fy s+2w i)£ 


(B36) 


Now the expansion (B35) is introduced under the integral of equation (B34) and the current 
II c (4, Wj) modified by the addition of the terms 2y£ + 3y 2 | 2 + . . . , is called hence- 
forth, I 2 c (£, w^. Thus 


I 2-C (5. w i) = e 


IWW . lT 2 


+ jx 2 f » r) ■ \ c ({, ’ w i ) d «’ 


+ 2 yjx 2 P r • k(j - ?)i 2 ;C (r,Wi) + 3 y 2 s • x 3 , . . (B 37 ) 


and K(£ - r) is given by equation (B17). Because y < 1 and equation (B37) is a linear 

integral equation, the last terms in equation (B37) may be expected to introduce only a 

small deviation to the solution (B36) already found, for I. . One may then assume 

1, u 

sI i(^ w i) + AI i(^ w i»y) 


j/3 a 0 w. £ o ri r- -7 

= e 0001 +jx 2 JT m- io[i 1 (r,w i ) + Ai 1 (r,w i ,y)l ar 


£ 

+ 2yjj? JT - r)[t l (€ , ,W i ) + dr (B38) 

Subtracting equation (B34) from equation (B38) gives 

AI 1 (?,w.,y) = A 2 K($ - r) * AI :l (r,w.,y) dr 


+ 2yx 2 Jr + AI 1 (r,w i ,yj] dr (B39) 
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In the last term of equation (39), the contribution of AI^ may be neglected under the 
integral as small compared with 1^ itself. Also, in equation (B39) (which is a correction 
in itself), a simplified expression may be used for the kernel from equation (B17); that 
is, for Oq© o 0 o « 1, K(£ - £» ) « j0 o /3 Q (^ - £')• After introducing 

\ < B40 > 

the following is obtained for equation (B39): 

ri A 

^l(^> w i ,y) = 2yA a J (?- !') Mid’, Wjyy) dr (B41) 


If the expression (B36) is taken for Ij (£, w^), is introduced into the first integral of 
equation (B34), and the first term is called <p(%, w.,y), one obtains: 


<ptt, w if y) = 2A a y 


J res - r)[ee 


■i w AV 


+ (i 


x jSCT 0' 

B)e 0 


j © r 

00 s 


i' 2 Wo w i« 


r 


dr (B42a) 


With an error of the order of 1 percent, B w 1/2 

-3 

tuted for the more important range of a Q < 10 
fication of equation (B42a) can be obtained: 


i 2a A w i 0 o* 

and e 0 ° 1 ° = 1 can be substi- 
and w^ < 10 and a significant simpli- 


, s 2yA | 
'" (| ’ w i )= 5 — j 

‘Vo 


/ »#»„! 


a cos a da 


2yA, 


a 





cos a da 


= - 2yij 


2 sin X^0 O | 

. *$>o« 


- 1 



(B42b) 


Note that terms with w^ dropped out of equation (B42b) because of the aforementioned 
simplifications . 
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The integral equation (B41) for Al^d, Wj,y) now assumes the form 


AIjU,y) = -2y| 


2 sin v g 


- 1 - cos vn + A, 


a 




^')Al 1 (r,y)dr (B41a) 


with v = xy/3 o 0 o for reasons of compactness. 

The integral equation (B41a) can be solved with the Laplace transformation method 
in analogy to the solution of equation (B15), preserving the same analytical approach to 
the entire problem. Thus, 


F v (p)= T Al^yjd* 


~\ 


2 ..2 


F„(p) = - *L - JL. - + 4 + 2y .. P— 1 - + V>) • 4 


P 2 + ^ P 2 * (P 2 + y 2 ) 4 


(B43) 


F„(p) = 2yv 2 jLzlL. 


(P 2 + 


Equation (B43) has two purely imaginary poles of the third order, p^ = + jy and Pg = -jv. 
Their residues are 


ResI >W eP5 - \ ^ ^ 

The Laplace inversion iT^[F i; (p)J leads to 

Aid, y) air ^^(p)] + 2 R es F^(p)e Pnl 

n=l 

2 

= sin i/f - y sin cos 

2v 2 2 


(B44) 


(B45) 
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The expression for Ig (£,w.,y) = 1^ c (|,w.) + Alj^ (eq. (B37)) introduces the effect 
of convergence, by AI^, on noise propagation in addition to the acceleration effects 
through a Langmuir -Blodgett potential field. This expression was described earlier in 
equation (B36) and is now written as 


I 2 (^ w i»y) = Be 


-^o°o 0 o^ 


+ (1 - B)e 


+j%i 3 o(s+ 2w i)0 o l 


+ AI(|,y) 


(B46) 


Since the imaginary part of equation (B46) is much smaller than its real part, the effect 
of convergence, by AI^, is to decrease the noise current Ig. For small values of I3 Q 
(i. e. , for a slightly nonlinear potential profile according to assumption (2) (p. 40) and 

. o Q 

< 1), the approximate expression AI is -(y/3)y | ; for small y, AI « 1. 
Thus, to a first approximation, this decrease is small. Incidentally, exact computations 
show that a convergent beam produces less noise current Ig than its nonconvergent 
counterpart, provided that both move through the same potential profile and that the non- 
convergent beam ends with the same X at the anode aperture. Note, however, that this 
conclusion is valid only for an electron beam of constant longitudinal temperature, as 
may be the case with electrostatically focused klystrons. 
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APPENDIX C 


NUMERICAL SOLUTION OF VOLTERRA EQUATIONS 
Statement of Problem 

Two Volterra-type equations are solved for two unknown, as well as two kernel, 
functions, all of which are complex functions of two real variables. The constant - 
temperature case CT and the variable-temperature case VT differ chiefly in their 
kernels. In both cases, the equations take the forms 

J(t w) = f(t w) - j K 1 (t, s)J(s, w) ds for J(o,w) = 1 (Cl) 

•'O 1 

U(t,w) = w • f(t,w) - J(t, w) + j f t K 2 (t, s)J(s,w) ds (C2) 

■•'0 

f (t, w) = exp|j© 0 J^(t, o) - 3(t,w)jj for 0 < t < 2 and w > 0 (C3) 


In the CT case, 0 is found from 


0 t (t,w) = ■ - 1 0 < t < 1 




f v(l) + 2ow 




Vl + 2ow 


1 <t < 2 


where t9(o,w) = 0 (see eq. (35)). 
The kernels have the forms 


K^t, s) = L(s)A r (s) 1 - Q r (t, s) 


K 2 (t,s) = L(s)A R (s)Q 2 (t,s) 


(C4) 


(C5) 

(C6) 


where 
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L(s) = 


and 



(1 - y) 2 R 2 (s) 


Qg(t, s) — 




exp(-w)p(t,w) 

f(s,w) 


dw 


See equations (30), (38), and (39). 

In the VT case, A R appears in the 3 equations but not in the kernel equations; that 
is, 


3 t (t, W ) = 


4 


0<t < 1 


”\ 


v(t) 


v(l) + 2crwA R (t) 
1 




1 <t < 2 


^1 + 2owA R (t) 

Kj(t, s) = L(s)[l - Q 1 (t,s)] 


J 


K 2 (t,s) = L(s)Q 2 (t, s) 


(C7) 


(C8) 

(C9) 


See equations (30) to (33). 

The problem is, in either case, to solve equations (Cl) and (C2) for J and U from 
t of 1 to 2 for w from 1 to 20. Then, use equations (CIO) to (C13) to solve for 

S(t) - n(t) : 




| J(t, w) | 2 exp(-w) dw 


(CIO) 


m-f 


|U(t,w)| 2 exp(-w) dw 


(Cll) 


n(t) + jA(t) = f J (t, w)U*(t, w)exp(- w) dw 
‘'O 


(C12) 
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S 2 (t) = - A 2 (t) 


(C13) 


Each case is either CT or VT. The value of a is 0. 001, 0. 0001, or 0. 00001; of y 
or a Q /r c , 1/2 or 2/3; of 0 Q , 4 or 8; and of X 2 , 0. 001, 0. 01, or 0. 03. The A R function 
greater than one uses m - f pairs of (2, 2), (2, 10), (3, 2), and (3, 10). 


Method of Solution 

Eight steps are involved in the process of solving the problem for a particular set of 
design values: (1) integrate the function 8; (2) compute the Q functions; (3) prepare 
kernels and 8 differences for the 0-to-l part of the problem; (4) integrate J and U 
from 0 to 1; (5) use these results to compute S(t) - n(t); (6) compute 8, the Q’s, the 
kernels, and the 8 differences for the l-to-2 part of the problem; (7) integrate J and 
U from 1 to 2; (8) use these results to compute S(t) - n(t). Each step consists of a 
separate pass on the digital computer. Steps (1) to (3) give card decks of output data; 
steps (4) to (8) have both input and output data in magnetic tape data files. 


Details of Numerical Technique 

The program uses two types of integration along curves in the complex plane. Where 
the chord, is close to the curve, trapezoidal integration is used, and where it is not close, 
spiral integration is used based on the formulas 



- sin A^) 



(C14) 



sin (a + bt) dt = -(cos Ag 


- cos A^) 



(Cl 5) 


where Aj =a+btj and j = 1,2. 

The program computes S - II as 1 - 211, even though direct calculation would 

2 

suffice for many cases of large cr and small X . That this is justified appears from 

o 

viewing a sample case with small ct (0. 00001) and large X (0. 01), where, at t = 2, 
II = -13.798, = 84 342, i// = 0. 79867, and A = -259. 17. This leads to ^$ = 67 362, 

A 2 = 67 167, and S 2 = 195. Now this conflicts with the condition S + n = 1, which is 
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known to hold everywhere, because if S is 1 - H or 14.798, its square is 218.98, not 

195. However, a slight change in 4>, and A (of less than 0. 01 percent) raises 
2 

by 12 and lowers A by 12, thus giving the correct S. The error in the calculation of 
these values is greater than this amount; hence, the altered mode of computing S - n is 
justified. 

The nature of the 3 equations near zero t for zero w or small w requires 
special care in the early integration. For this reason, 20 of the 70 values of t from 0 
to 1 are below 0. 02. The w values need to be finely spaced, near zero because of the 
nature of the -3 function and near one because of the exp(-w) weighting of certain 
integrals. The Q calculations use 101 values of w computed from 

w(n) = 0. 002n^ where n = 0, 1, . . . , 100 

The J and U integrations use only 48 of these values (chiefly to lower the time of 
solution), and they correspond to the following 48 values of n: 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

15 

17 

18 

20 

21 

22 

23 

24 

25 

27 

28 

30 

32 

33 

35 

37 

38 

40 

42 

43 

45 

47 

48 

50 

52 

53 

55 

57 

58 

60 

70 

80 

90 

100 


Note that these values are finely spaced near aw of 1, or an n of 22 and 23. 
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APPENDIX D 


POWER SERIES FOR LANGMUIR-BLODGETT POTENTIAL 

From reference 13, the following relation between V(x) and V(a) may be 
obtained: 


V(x) = V(a Q ) 




r 


«V»/V 


The a function was calculated by Langmuir and Blodgett in a series 


a = y . J-y 2 + Jly 3 _ 0. 1432 y 4 + 0. 00216 y 6 - + . . . 
10 100 


(Dl) 


In the present investigation, 


y 





In 1-| 


(D2) 


Since y = (a../r ) < 1, y itself may be expanded into a series and substituted in equa- 

U v 

tion (Dl). The result is 


^ = a^-2. ^ s a(y£) = -yt(l + |y| + lly 2 | 2 + 0. 662 y 3 | 3 

+ 0. 611 y 4 £ 4 + 0. 595 y 5 | 5 + 0. 57 y 6 | 6 + . . (D3) 

In equation (D3), the first three terms are exact; the others are rounded off to the nearest 
third digit. 
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APPENDIX E 


A 

A R 

a 

a o 

b 

c 

E 

e 

F,f 

^min 

G(w.) 

g(w) 

H(x, c) 

h(p) 

m 

i 

J(|) 

K(|, ?) 
k 

^ _1 |F(p)] 

M e 

m 

P M 

P 

R(x) 
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SYMBOLS 

area 

total area compression ratio 
shield radius 

distance between cathode and aperture 
beam radius 

|“a/( a o/ r o3^ or constants 

electric field 

IQ 

electron charge, 1. 6x10“ C 
distribution functions 
minimum noise figure 
integration constant, eq. (16) 
direct-current distribution function 
inhomogeneous part of eq, (11) 
complex function 

complex radiofrequency -noise convection current 
current density 

modified complex noise current, I(|)/A R 

kernel of integral equation 

Boltzmann constant, 1. 38xl0” 23 J/K 

inverse Laplace transform, si(£) 

Laplace transform, sF(p) 

31 

electron mass, 9. 11x10 kg 

bjb n 

a' o 

beam perveance, A/V 3 / 2 
Laplace variable 

plasma frequency -reduction factor 



radius 

unit step function 


spectral function of noisiness (normalized), |i//(§)$(£) - A 2 (£)j 


temperature, K 


complex electrokinetic noise voltage 

velocity 

potential 

dummy variable 

energy variable 

distance 

a /r 
o' c 

Langmuir -Blodgett function 

M e /2kT c 

M e /2kT c x A R = o 1 q /A r 
transit time parameter (appendix B) 
cathode opening angle 
In r/r Q (only in appendix D) 

Dirac delta function 

dielectric vacuum constant, 8. 854x10" 12 F/m = C/(J)(m) 

e/M e = 1.759X10 11 C/kg 

transit time parameter, wa o /( 27 jV a )^ 2 

transit time function 

beam size polynomial 

-A0 2 0 

noise power spectral function (imaginary part of) 

"p,a /M 

+x© 0 '\^ 

normalized distance, x/a 
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n(£) noise power spectral function (real part of) 
p space charge density, C/ cu m 

cr T kT c /2eV a xA R = a o xA R 

HZ) J fU(|,w i )| 2 e 1 dw. 

HZ) |l(^,w i )| 2 e 1 dWj 

oj frequency x 2n, Hz 

Subscripts: 

a anode or aperture 

c cathode 

e electron 

i injection class 

l longitudinal 

o refers to direct current or ultimate beam size 

p plasma parameter 

q reduced plasma parameter 

s slowest 

t transverse 

1, 2 radiofrequency quantities, first and second approximations 

Superscripts: 

~~ average 

( ) average over distribution 
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